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Abstract 

The behavior of a stationary inverted point mass pendulum pivoted at its lower end in a gravi- 
tational potential is studied under the influence of statistical fluctuations. It is shown using purely 
classical equations that the pendulum eventually tips over i.e evolves out of its initial position of 
unstable equilibrium, and, in a finite amount of time points down assuming a position of stable 
equilibrium. This 'tipping time' is calculated by solving the appropriate Fokker- Planck equation 
in the overdamped limit. It is also shown that the asymptotic time solution for probability cor- 
responds to the Boltzmann distribution, as expected for a system in stable equilibrium, and that 
the tipping time tends to infinity as the parameter corresponding to the strength of thermal fluc- 
tuations is tuned to zero, thereby defining the limit where one recovers the classical result that a 
stationary inverted point mass pendulum never tips over. The paper provides a unique perspective 
showing that phenomena like tipping that have been often attributed to quantum mechanics can 
be studied even in the domain of purely classical physics. 
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I. INTRODUCTION 



Inverted stationary pendulums have always been an interesting study. Classically, in 
the absence of any perturbations to the gravitational potential, such a position of unstable 
equilibrium corresponds to a saddle node in phase spaced If the system is conservative i.e 
there is no damping, these saddle nodes are connected by heteroclinic trajectories of infinite 
time period. In the presence of damping, these trajectories change to spirals and fall into 
stable fixed points. In both cases, however, as long as the system is left undisturbed, it 
stays precariously at the saddle node. Simply put, from such a classical standpoint, an 
unperturbed stationary inverted point mass pendulum stays inverted forever. The scenario 
changes when one seeks a quantum mechanical solution. The uncertainty principle forbids 
fixing both the position and velocity simultaneously. One can at best localize the initial 
position of the pendulum within a 'small enough' region around the inverted position and 
study its behavior. Work has been done in the past£~- to explore how the uncertainty 
principle causes quantum rods and pendulums to get delocalized and tip over. 

So far, this tipping has been well studied as a quantum mechanical phenomenon. Con- 
sidering a pendulum with friction, say in a jar of olive oil, it is hard to believe that tipping 
can happen only through quantum effects. Even in the classical world, the fluctuation- dis- 
sipation theorem tells us that in a dissipative system, there are fluctuations of a statistical 
rather than quantum origin which can very well push the pendulum out of equilibrium. 
In this paper aimed at undergraduate students, graduate students and teachers of physics, 
we return to the classical pendulum, study its evolution in the presence of finite tempera- 
ture statistical fluctuations and show that even classical effects can lead to such interesting 
behavior. The pendulum is classical in the sense that its initial position can be exactly 
fixed to the inverted position of unstable equilibrium with zero momentum. We work in the 
high friction (overdamped) limit and use the corresponding Fokker- Planck equation for this 
system to study its behavior and quantify the time it takes to tip over. 

II. FORMULATION OF THE PROBLEM 

Consider a point body of mass 'm' fixed to one end of a massless rigid rod of length '1', 
while the other end is pivoted. Let this system be free to move in a vertical plane under the 
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influence of a gravitational potential of the form: 



V{6) = mgl cos 9, 9 E [0,2tt] (1) 

where 9 is measured from the vertical. Now, we introduce white noise from statistical 
fluctuations in the form of a stochastic variable, £(t). The noise is white in the sense that it 
is delta correlated in time: 

(am(t'))=^k B T5(t-t') (2) 

where, 7 is the coefficient of friction, k B is the Boltzmann constant and T is the abso- 
lute temperature^. With these, we can write the 'Langevin Equations' of motion for the 
pendulum as: 

(me) m^ 10(t) _'im +m (3a) 

We work in the high friction limit and linearize the above equations by dropping the intertia 
term to get, 

Using Eqs. HI [I] and [2], we can expect to obtain a formal solution to the expectation value 
(*(*)>■ 

Equivalently, we can work with probabilities. Defining p(9, t) as the probability density 
of finding the pendulum at position '#' at time and using Eqs. H]and|2]in the continuity 
equation associated with p(9, t), we can find the evolution of the expectation value (p(9, t)) = 
P(9, t). The 'Fokker- Planck' equation (corresponding to Eq. Hj) that describes this evolution 
can be written as: 



dP(9,t) 1 d (dV{9) m dP{9,t) 



(5) 



The steps involved in moving from Eq. H]to Eq. [5] can be found in Ref |5|. We prefer to work 
in the probability picture and perform our analysis with Eq. [5j 
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III. NUMERICAL SOLUTION TO THE FOKKER- PLANCK EQUATION AND 
ESTIMATION OF THE TIPPING TIME 

Substituting Eq. [T] and denning the dimensionless quantities r = — —t and a = 3 we 

7 mgl 

can rewrite Eq. as: 



dP(0,r) 8P(9,t) . . d 2 P(9 



T 



()T P(0,T)cceO-—^Bine + a—^. (6) 

Henceforth, we shall work with the dimensionless time V and any reference to 'time' would 
mean the latter. Also, 'a' can be regarded as the strength of thermal fluctuations. We 
assume that at r = 0, the pendulum is exactly at 9 = i.e: 

P(0,0)=6(0) (7) 

Also, since the pendulum moves in a circular space, we impose periodicity in 9: 

P(9 1 r) = P(9 + 2n 1 r) (8) 

Using Eq. El with the boundary conditions Eq. [7J and Eq. El we can obtain a solution for 
P(9,r). 

To do this analytically we would have to obtain the solution as an eigenfunction expan- 
sion by spectral resolution of the differential operator on the right hand side of Eq. |6j First 
however, the operator has to be made self adjoint by using a suitable transformation on the 
independent variable. A detailed prescription to do this can be found in Ref. |5|. Such a cal- 
culation is mathematically cumbersome. We prefer to simply solve the problem numerically 
and find that it captures the essence of the solution effectively. 



A. Asymptotic time solution 

We can calculate the solution (P(9, <x>) = $(#)) for Eq. [6] which is the probability dis- 
tribution when the system has completely relaxed to its stable configuration, by setting the 
time derivative to zero: 

„ d&(6) „ d 2 <$>(9) 
-$(#)cosfl-— - ^sinfl + a — ^ = (9) 
do at) 

As one expects for a system in stable equilibrium, the solution for Eq. considering Eq. [7] 
and Eq. [8] is a normalized Boltzmann distribution centred around 9 = n. 



where Io is the zero order modified Bessel function of the first kind and the prefactor comes 
from normalizating <&(9) on the circle. We expect that long after tipping, the probability 
distribution assumes the form of $(#) displayed in Eq. [TUJ 

B. Evolution of probability density 

Fig. [H shows the evolution of the spatial probability distribution with time. It can be 
seen that with time, probability oozes out of the initial delta peaks, the maximal bumps 
travel towards 9 = n and upon reaching it, stay there. It is also seen that the asymptotic 
probability density indeed corresponds to of Eq. [10J Another noteworthy observation 
is the symmetry of the solutions about 9 = tt. This makes physical sense as the pendulum is 
equally likely to tip clockwise as it is to tip anti- clockwise and during the actual process of 
tipping, the pendulum breaks this symmetry. Mathematically, the fact that the probability 
solutions respect parity can be confirmed by observing that Eq. [6] is invariant under the 
transformation 9 — > —9. 

Fig. [2] shows the time evolution of P(9, r) at various angles. It is seen that the graphs do 
not oscillate which is typical of solutions in the high friction limit when the inertia term is 
neglected. Once again, we see that at large times, the probability density is maximum for 
9 = 7i. This confirms that after a very long time, the pendulum is most likely to be found 
at the stable equilibrium point of 9 = tt. One can thus define the tipping time, r tip as the 
time when the probability is first a global maximum at 9 = n. 

C. Calculating the tipping time and recovering the classical mechanical limit 

Defining the tipping time as the instant when probability is first maximum at 9 = n 
suggests that the value of T tip can be calculated by solving the equation 



However, observing Fig. [T] we see that 9 = n is always a point of local extremum and 
therefore such a slope is always zero. This arises from the symmetry about 9 = n mentioned 
earlier. On the other hand, we can see from Fig. [T] that as 9 = ix makes a transition from 
being a point a local minimum to global maximum, the second derivative changes sign from 








(11) 
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positive to negative. Thus, we can recognize Tu p as the instant when the second derivative 
becomes zero, by solving the equation 

Fig. [3] shows the variation of r tip with a. It can be seen that when a — > oo i.e when the 
thermal fluctuations greatly dominate, the pendulum tips instantaneously (tu p — > 0). On 
the other hand, when a — > 0, we find that Tu p — > oo which is tantamount to saying that 
in the limit of negligible thermal fluctuations, we recover the limit in which the pendulum 
never tips over. 

Finally, let us calculate the tipping time for a typical system to get a feel for the order 
of numbers involved. Considering a point mass pendulum of length I = 10/xm with a 
spherical pollen grain of mass m = 10~ 14 kg and radius r = l/j,m acting as the point mass, 
immersed in olive oil with coeffecient of dynamic viscocity rj = O.lPa.s at room temperature 
(T = 300K), we can calculate T Up for this configuration which corresponds to a = 0.0042 as 
3.206 dimensionless units. Comparing the definition of 7 in Eq. [Hand the formal definition 7 

of 77, it is easy to see that the two are linked by 7 = r](nr 2 )l. With this, we can calculate 
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the observable tipping time t tip = -Tu p as 21.919 seconds. 

mgl 

IV. SUMMARY AND CONCLUSION 

The evolution of an initially stationary and inverted classical pendulum in the presence of 
statistical fluctuations is studied and it is found that these fluctuations push the pendulum 
out of unstable equilibrium. The time taken for the pendulum to tip to the configuration of 
stable equilibrium is calculated numerically and it is seen that within the confines of classical 
analysis, only in the limit of negligible thermal fluctuations does the pendulum stay inverted 
forever. Analysis of this simple system elegantly demonstrates physical phenomena that 
commonly occur in dynamical systems out of stable equilibrium and hence can be used not 
only in pedagogy but also as an effective toy model 6 to study complicated systems and their 
behavior, in wider and more profound contexts. Also, the analysis being purely classical 
demonstrates the abundance of highly interesting physics in the relatively forgotten realm 
of classical physics that are still unexplored and are worthwhile to study to gain a rich insight 
into the physical laws governing observable reality. 
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Figures 




FIG. 1: Evolution of P(9,t) versus 9 distribution setting a = 0.5 



P(0 ,T) 




FIG. 2: Evolution of P(6, r) at various setting a = 0.5 
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FIG. 3: Tipping time, r t i p versus thermal strength, a 
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